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a SUMMARY 
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| The present work is a direct continuation of an earlier investigation by the author on the 
e subject [1]. The question of finiteness for the energy per unit cell of the crystal has been 
ussed. A numerical mistake, which, however, does not change the main qualitative con- 
isions made in [1], has been corrected. The energy bands have been computed for five differ- 
nt choices of the d orbital exponent. The result of thes2 new calculations indicates that a 
markable increase in the d orbital exponent, when chosen according to the Slater rules, is 
squired to minimize the energy. 


Introduction 


by Kiessling [2] and theoretically by the present author [1, 3, 4] and by Longuet- 
gzins and de Roberts [5]. In the present investigation the secular equations 
or the different symmetry adapted orbitals of the B, valence electrons are solved 
1 the spd approximation with variable d orbital exponent. 

Section 1 gives the effective Hamiltonian operator for an electron in the average 
d of the rest of the crystal. The total Hamiltonian operator of the system 
sidered has been given and divided up into equal parts comparable to the 
tive operator. It turns out that both operators are finite and have the same 
er of magnitude. Section 2 gives the numerical results of the present calcula- 
ns including variations of the d orbital exponent. The existence of conduction 
mds and the great influence of M shell orbitals strongly dependent on their 


The effective Hamiltonian operator of CaB,1 


We denote the valence electrons by 7 or j, the boron nulei (nuclei + inner 
rons) by B, or B and the ionized Ca atoms (considered as point charges) 


1 Before discussing the Hamiltonian operator it may be convenient to correct some misprints 
ing in [1]: 


- (By) 
530, eq. (37). A factor 2 is missing for the term = 1/r;- (The notation By here stands for the 
s “belonging to” the electron 7). 
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by M, or M (metal). The total Hamiltonian operator for the system under con- 
sideration is given by eq. (37) of [1]: 


Zu Zu 


M, M+M, ym, BM 1pm 


Le ig if (1) 


© 


The constant nuclear energy in the second line of (1) must be added to get a 
finite crystal energy per unit volume. This also implies that the cohesive energy 
with opposite sign is given by (1). The nuclear charges are 


Zn=Zp,=3; Zu= Zu, = 2A, (2) 


where A=1 for complete ionization of Ca. 

Our intention is now to split the Hamiltonian operator (1) into an infinite 
sum of identical formal one-electron operators. We denote by hy the Hamiltonian 
operator of one electron in the field of all the nuclei: 


ENDS tS 


Brg M ly 


1 


The potential energy part of (3) tends towards infinity for an infinite crystal, 
since the number of atom neighbours at a certain distance is proportional to 
the squared atomic distance. 

The coulomb interaction energy between one electron and all the others may 
be written 


gee (4) 


i+111i 


The nuclear B-B interaction energy for one boron atom divided by the average 
number of valence electrons per boron atom, 


can be written ; 
Zs 1 3 
res “8 
2 Zp+Zy/6 B+B, TRB, 1 +4/9 2. (6) 


Here we have introduced the notation 


(8) 
p- 531, eq. (38). A factor 4 is missing for the term ¥ 2/r wp. 
M 0 


Danse: eq. (48). The term, Z,/R, must be replaced by Z,/R,. The sign of the term containing 
x must be changed. 
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This term is also infinite. Similarly, we can introduce the nuclear M-M inter- 
action energy “‘per valence electron” 


By 1 243 1 Aa I 
+B: (8) 
6Z,+ Zu ™M Tym, 9+AM rym, 9+163 


Eyes 


where f is a finite positive constant with a definite value for 4=1 (8 ~ 0.2 a.u.). 
The electron-electron energy (4) can also be expressed by use of (7): 


he=(3+ 4/3) 3 +6, (9) 


where 6 is a finite operator with expectation values determined essentially by 
the electronic interaction between one valence electron and all the other valence 
electrons within the same boron atom (order of magnitude 6~0.6 a.u.). 

Furthermore we introduce the average nuclear B-M interaction energy per 
electron, 


Z; Zu 1 2A 1 


= = ‘ 10 
oS Z3+Zy/6 mu 1am 1+A/9 MY?pmu ( ) 
Finally the ionization energy, H;, may be taken per electron, yielding 
ve 
= SS 1 
PRETTY (11) 


where J is the ionization energy of one Ca atom. 
The total Hamiltonian operator (1) ‘‘per electron’ may then be written by 
using (3), (4), (6), (8), (10) and (11) as a formal “one-electron operator”: 


hy=ho tthe th est teu t Emu t €. (12) 


Since, however, hf, in fact is a two-electron operator, (12) cannot be used as an 
effective one-electron Hamiltonian operator. To construct such an effective opera- 
tor, we must compensate for the electron-electron interaction by using effective 
charges, Z, chosen differently for different atomic distances and different atomic 
orbitals involved in the integrals occurring for the energy expectation values 
(see [1j, Table 6). | 

The effective Hamiltonian operator used for the energy band calculations in 
[1] is given by [1], (38): 


ff " eff eff 9 
hes EA ae eee, 1 ( 23 5 Bay )+en 13) 
M 


a at - > 
2 Srp My 1+A/9\Bren, M Tap, TMB, 


The first three terms of (13) are dependent on the atomic orbitals involved. 
Here, Z# is chosen to satisfy the relation 
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—-A-¥—-5— why th, (14) 
a M ‘yy 


Baty 


as well as possible. The nuclear energy in (13) added to this quantity in order 
to get a finite energy, is the interaction energy between one boron nucleus and 
all the other atoms per electron. Since the valence electrons shield the boron 
nuclei, effective charges occur also here. These effective charges tend indefinitely 
towards Z°*= —4/3 for large distances. Consequently we have a relation similar 
to (6): 


1 w2¢ 1/3 
eff __ 29 — + oy, 15 
SF WEP MOET: 1+1/94 x (15) 


where « has a small finite value (c=0.02 a.u. for 2=1). 
Using the equations (14), (15) and (10) we can write (J3): 


h~hothe+ e% + eau + 1. (16) 


Since the effective charges tend towards —A/3, the sums for large distances can 
be calculated as pure Madelung energies. As is well known, these energies are 
convergent under the condition that the terms are added in suitable sequence 
of neutrally charged groups [6]. 

Let us now consider the difference between (16) and (12): 


h—h,~ thet ef — 3 ep —} eq. (17) 


Inserting into (17) the relations (9), (15), (6) and (8), we find that the coeffi- 
cient of > vanishes identically, giving a finite difference 


b& 


h-h,w~at4(6—B)~0.2 au. (18) 


Thus both h, and h are finite and their expectation values have the same order 
of magnitude. 

The effective Hamiltonian operator (13) has been used for determining the energy 
bands. The sums over neighbouring atoms have been performed for the 24 nearest 
boron atoms and the 8 nearest metal atoms. The remaining terms have been 
summed up as electrostatic point charge interaction (Madelung energy) and put 
in the last term of [1], (38). The notation e, is used for this quantity per B, 
octahedron. : 

The quantities (47) and (48) of [1] entering in the secular equations, which 
give the energy bands corresponding to (13), can be written (with the approxi- 
mations and notations introduced in [1)): 


(24) A, (8) | 
ho =< By| h| Bo> =Co— Zp By > Zp (p, - 3") —> 2A (p,- 22) +€5 Ag; (19) 
b m m 
on Zp 
h=C,-Z, By ot feats (20) 
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Ey (14, Ey(it,,) 
0 kk = (0,0,0) Ik = (0,",m) 


mt 1 
K = (oygr2 


Fig. 1. Energy eigenvalues, éx, given in a.u. for the bonding orbitals for some different wave 
vectors in the LM approximation as functions of the d orbital exponent Z/3. The straight lines 
correspond to the L approximation. 


where ¢) is a constant not depending on the special orbitals involved: 


a % (24) Zs ' 3) 
- —2 = — I}. : 
Eq seit? ER 15 Fs 35 (21) 
For complete ionization (A=1), e,= —0.059 a.u. The negative sign of this quan- 


tity indicates that the stability of the crystal is somewhat increased on account 
of ionization. 

Jn determining the one-electron crystal energy we have not compared it to the 
mean binding energy of an electron in a free boron atom, which, however, might 
be more realistic than to compare it to a state where all the considered par- 
ticles are free. Thus this binding energy term should also be added, increasing 
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Fig. 2a. Energy bands é¢ in the L approximation. k= (9, 0, 0). 


the orbital energies somewhat. This might be taken into account by increasing | 
the constant term é . 

The energy bands given in section 2 have been calculated, for simplicity, by 
putting e,=0. To get the actual bands corresponding to (13), the given bands 
thus ought merely to be lowered by the quantity 0.059 a.u. But, since the 
zero point level cannot be meaningfully defined with very great accuracy, this 
common lowering of the bands has no certain importance and can be neglected. 

Unfortunately, the mistake of sign in [1, (48)] also appeared in the numerical 
calculations, but renewed calculations show that the main qualitative conclusions 
still hold. In addition to these calculations variation has been performed for the 
d orbital exponents, the result of which will be given in section 2. 
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Q.u. 
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Fig. 2b. Energy bands && in the L approximation. k= (0, y, 0) and (0, 0, 7). 


2. Electron distribution and energy bands for different d orbital exponents 


The molecular integrals occurring in (19) and (20) have been calculated by 
means of a programme, made by the present author, for the electronic computer 
Besk [7]. This programme is free for anyone to use and a special table com- 
puted by use of it can be’ bought [8]. Also the different secular equations have 
been solved by means of Besk. 

The electron distribution and the energy bands were given in [1] in the D 
approximation (including hybrids of 2s and 2>p orbitals) and in the LM approxi- 
mation (including hybrids of 2s, 2p, 3s, 3p and 3d orbitals). The Z shell orbi- 
tal exponents were chosen according to Duncanson and Coulson [9]. For the 
M shell, the orbital exponents were chosen according to Slater’s rule, yielding 


the effective charge Z=1. 


at 
or 
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0 n/2 ov 


Fig. 2c. Energy bands é; in the L approximation. k= (9, , 0) and (9, 0, (7). 


The overlap between M shell orbitals of neighbouring atoms is, however, very 
great, indicating that the orbital exponents should rather be chosen according 
to the conditions in the whole crystal than in the free atom. This great overlap 
also causes great overlap elements for the corresponding crystal orbitals occur- 
ring in the secular equations. Thus the eigenvalues and eigenvectors are very 
sensitive to small changes in the secular elements. Moreover, very many neigh- 
bours must be taken into consideration on account of the long range interaction 
for M shell orbitals with small orbital exponents. Since it would lead to quite 
extensive calculations to take sufficiently many neighbours into account, some 
of the secular overlap matrices are indefinite in the present approximation. 

An increase of the effective charge of the M shell orbitals to about Z=5 would, 
however, decrease the interaction range for these orbitals to the same order of 
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-0.5 


Itoy 
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Fig. 2d. Energy bands éy in the ZL approximation. k= (0, 9, ~). 


magnitude as for the LZ shell orbitals (cf. [7], Figs. 2, 3 and 4). The secular 
equations corresponding to the interrupted series of neighbours in the atomic 
integrals (cf. (19)) has consequently in general no indefinite overlap matrices. 
For somewhat smaller choice of Z, however, some of the overlap matrices are 
indefinite, or nearly so, and the corresponding secular equations cannot be solved 
without taking many more neighbours into account. 

Since 2s and 3s orbitals are only slightly different, when the orbital exponents 
are chosen to give nearly equal interaction range, t.e. the overlap between nor- 
malized crystal orbitals constructed from 2s and 38 orbitals tends towards 1, 
corresponding to a ‘‘nearly indefinite’? metric, it gives very little extra informa- 
tion to make use of both. The same thing is valid for parallel 2p and 3p orbi- 
tals. Thus 3s and 3p orbitals have been excluded from the present calculations 
and the M shell has been permitted to contribute with d orbitals only, the 
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Fig. 2e. Energy bands éy in the L approximation. k= (9, 9, (). 


symmetry of which is very important on account of the point group symmetry | 
of the crystal. 

Five different choices for the effective charge, Z=1, 2, 3, 4 and 5 have been 
made for the d orbitals in the present investigation. The modification of the d 
orbital exponents makes it necessary to use also modified effective charges for 
those integrals of type (19) and (20) of the effective Hamiltonian operator, con- 
taining d orbitals. The calculation of these new effective charges has been made 
according to the equations (49)-(53) occurring in [1]. Thus [1], Table 6 can 
be widened for the effective charges Z) of the diagonal one-centre energy inte- 


grals containing 3d orbitals. For the case 4=1, the result is given by Table 1 
for different d orbital exponents Z/3 (Z=1, ... 5). 
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Fig. 3. Regions of the energy bands é« for the bonding and conducting orbitals in the Z and 
IM approximations. Dotted lines correspond to rather uncertain LM bands not occurring in 
the L approximation. 


Table 1. Effective charge, Z,, for the diagonal one-centre integrals as a function of 
the d orbital exponent Z/3. 


4 
0.9697 


5 
1.1317 


2 
0.7158 


3 
0.8213 


Z 
Zo 


1 
0.6700 


As we see from Table 1, Z, for 3d orbitals with Z=5 is closest to the corre- 
sponding values for 2s and 2 p (1.3540 and 1.2607, resp.), which is quite reasonable. 

Considering the Z dependence for diagonal 3d integrals of (19), which are in- 
dependent of the different types of 3d orbitals, we find a descending curve, mainly 
on account of the increase of Z). The most essential contribution, however, comes 
from the two-centre terms (20). It is found that the secular elements are rather 
considerably affected by the variation of Z. Certainly, the variation interval 
should be chosen much narrower to make it possible to draw reliable diagrams, 
but since a very great part of the calculations must be performed by use of 
desk machines, the time required to make these calculations would be quite long. 
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Table 2. Energy eigenvalues of the effective Hamiltonian operator (13) in a.u. for tha 
occupied symmetry adapted orbitals in the L approximation (first line) and in the 
LM approximation (the following lines) given for different points in the wave vector 
space. On account of some indefinite overlap matrices several values have beer 
omitted or put within parentheses. 


Approx. | 


L 
IM Z= 


oF WNW Re 


LM Z= 


oO FR whe 


LM Z= 


Coe & bo = 


LM Z= 


oF WON RE OR wD 


IM Z= 


oR whe 
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— 0.916 | — 0.850 
— 1.072 | —1.143 
(— 1.576) 
—1.869| — 1.514 
—1.211] —1.215 
— 0.939} — 0.992 
(— 0.956) 
(— 2.330)|(— 1.656) 
— 1.370] — 1.583 
— 0.566 | — 0.572 
— 0.646 | — 0.627 
— 0.714] —0.714 
— 0.659 | — 0.642 
— 0.592 
— 0.575 | — 0.578 
— 0.166 | — 0.137 
— 0.167] — 0.140 


(— 0.280) (— 0.222) 
(—0.296)|(— 0.185) 


— 0.220] — 0.148 
— 0.179] — 0.140 
0.058| 0.066 
(0.055)| 0.065 
(0.012)](— 0.110) 

— 0.023 
0.031] 0.040 
0.268] (0.221) 


— 0.923 
— 2.762 


— 2.955 
— 1.743 


OO) 


(— 1.274) 
— 2.026 


— 0.583 
— 0.607 
— 0.820 
— 0.628 
— 0.592 
— 0.585 
— 0.086 
— 0.091 


— 0.088 
— 0.086 


0.075 
(0.075) 


(— 0.107) 
— 0.009 
0.051 


(0.113) 


(— 0.341)}(— 0.120) (— 0.068) 
(— 0.222),(— 0.129) (— 0.079) 


3.176 
0.492 


(0.500) 
(0.584) 


0.477 


(2.489)| 
3.539 


0.629 
(0.558) 
(0.635) 
(0.702) 


0.626 


(2.030) 
4.438 


0.816 
(0.758) 
(0.820) 
(0.875) 


0.814 


— 0.957 
— 0.980 


— 1.874 
— 1.316 


— 0.917 


— 2.105 
— 1.227 


— 0.572 
— 0.627 
— 0.714 
— 0.642 
— 0.592 
— 0.578 


— 0.141 
— 0.141 


— 0.185 
— 0.151 


— 0.015 


— 0.170 
— 0.073 
— 0.042 


0.154 
— 0.238 
0.222 
0.737 
1.618 


0.128 
0.115 


0.029 
0.066 
0.104 


— 1.017 


(— 2.017) 
— 1.972 


— 0.899 
— 1.681 
(— 2.466) 


— 1.345 
— 1.132 


— 0.583 
— 0.607 
— 0.820 
— 0.628 
— 0.592 
— 0.585 


— 0.102 


— 0.146 
= 0.199 
— 0.138 


— 0.061 
(— 0.068) 


(— 0.067) 
— 0.094 
— 0.080 


0.075 
0.041 
0.236 
0.546 
1.198 


— 0.002 
— 0.003 
— 0.040 
— 0.086 
— 0.101 
— 0.046 


7 
io > 0, { 
ova 
— 0.878 | — 1.009} — 1.014 
— 0.991 
(— 3.939) 
— 1.727 |(— 3.814)} — 1.940 
— 2.803 | — 1.470 
— 0.965 | — 1.003 | — 0.902 
— ead 
— 1.647} — 1.892 
— 1.766 | — 1.456 
— 1.370} — 1.412 
— 0.582 | — 0.628 | — 0.582 
— 0.621 |(— 0.628) — 0.621 
(— 0.657) 
— 0.626 | — 0.629 | — 0.626 
— 0.592 — 0.592 
— 0.585 | — 0.629 | — 0.585 
— 0.109 | — 0.012 | — 0.107 
— 0.111 |(— 0.024) 
(— 0.104) 
(— 0.230) 
— 0.119} — 0.034] — 0.144 
— 0.016} —0.115 
— 0.010} — 0.053 | — 0.073 
(— 0.071) (— 0.076) 
— 0.266 (— 0.256) 
— 0.069 | — 0.093 
— 0.036 | — 0.073 | — 0.094 
(0.905)| (0.288)} 0.120 
(—0.109)} _(0.147)) — 0.049 
(0.887)} 0.246 
0.936 0.737 0.597 
1.771 1.378 1.267 
0.169 0.042 0.074 
(0.159)} (0.040)} 0.073 
(0.144)} (0.019) 
0.157 0.003 | — 0.078 
0.130 | — 0.028 | — 0.047 
0.153 0.010 0.030 


er eee | 
; a L 1 m1 7 7G Te 
Wave vector — | (0, 0, 0) G 0, 0) (xt, 0, 0) (0. > 0) (0, , 0) (=. 7 0) (xt, zt, 0) (0, ry °) (0, 7, 2) (5, Port 


7 at 
(0, 0, °) (0, 0, 7) (5, 0, 
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Table 2 (continued). 


ector > | (0, 0, 0) ( 0 0) (xt, 0, 0) (0. a 0) (0, 2, 0) (= = 0) (x, 2, 0) (0, ide eS (5, a >) ee 
Approx. (0 0 ) (0, 0, qt) (= 0, 5) (2, 0 qt) 
Ty, 0.082 0.017 | — 0.032 0.108 0.143 0.031 | — 0.016 0.164 0.452 0.064 0.075 
z= 1 0.062 0.013 | — 0.034 0.089 0.028 |(— 0.023) 0.147 (0.413) 0.060 (0.059) 
2 (0.042) (0.431) 
3 |(— 0.051)} — 0.052 — 0.071 | —0.137 (0.536)| — 0.095 |(— 0.126) 
4} 0.030} — 0.031] — 0.091} 0.058] 0.083] —0.019| — 0.085] 0.115] (0.442)} 0.009] —0.049 
— 5} +0.073} 0.008] — 0.046] 0,096} 0.125} 0.019] —0.038] 0.148] 0.441| 0.049] 0.033 
M Z=1 (0.203) 0.673 (0.304)} (0.442) 3.074 (0.355)| (0.966)| — 0.017 0.163 
2 (— 0.149) (— 0.183) (— 0.181) (— 0.096) (— 0.376) — 0.092 — 0.088 
3 |(— 0.462) ;(— 0.601) (— 0.362)| — 0.021 (0.005)| — 0.086 (0.122) 0.050 | — 0.008 | — 0.009 
4 0.431 0.634 0.548 0.309 0.371 0.452 0.665 
5 0.864 1.028 1.295 0.818 0.793 0.945 1.245 0.773 0.724 0.872 0.903 
oa 
Table 3. Coefficients occurring in the crystal orbitals of the occupied energy bands 
g MA Pp g) 
for the atomic orbital parts of the octahedral site orbitals given by Fig. 4 as func- 
tions of the wave vector. The coefficients are given in the JZ approximation and in 
if the LM approximation for Z = 3 for a few wave vectors. 
(0, 2, 0) (x, 2, 0) 
i? (0, 0, 0) (a, 0, 0) (0, 0, 2) (20, 0, 7) (0, 2, 7) (7, 7, 7) 
a a 
nation} L CA 4 aE Le LM L LM i LM : LM 
- - 
vital 
EF 0.981 0.932) 0.925 0.847} 0.976 0.912 0.908 0.930 0.965 0.936 0.880 0.827 
i 0.109 — 0.030} — 0.273 —0.219] 0.096 0.097| — 0.296 0.014 0.065 0.120] — 0.327 0.138 
4 — 0.161 0.034] — 0.263 — 0.002) — 0.195 0.129] —0.297 —0.355| —0.254 —0.003|] — 0.344 0.395 
ido 0.325 0.343 0.328 0.087 = 0.215 — 0.028 
da 0.153 0.342 0.187 0.049 0.251 0.375 
is 0.933 0.694) 0.921 0.977| 0.949 0.961 0.939 0.957 1.000 0.995 1.000 0.996 
po |—0.361 0.413] —0.388 —0.038]—0.315 — 0.248} —0.343 —0.113] —0.011 — 0.042] — 0.020 0.059 
do — 0.561 ~ 0.033 0.041 ~ 0.245 0.071 0.013 
— 0.180 — 0.209 — 0.118 — 0.104 — 0.054 — 0.066 
0.989 0.916} 0.974 0.948} 0.974 0.948 0.948 0.954 0.948 0.954 0.893 0.853 
0.146 —0.160} 0.227 —0.012] 0.227 — 0.012 0.320 0.295 0.320 0.295 0.450 0.497 
I — 0.367 — 0.317 — 0.317 ~ 0.044 —~ 0.044 0.160 
> |—0.405 —0.209|—0.693 —0.468|—0.407 —0.315] —0.482 0.850] — 0.691 —0.168] — 0.601 —0.307 
>a | 0.532 0.479] 0.333 —0.306] 0.531 0.410) 0.460 —0.301] 0.254 0.438] 0.339 — 0.176 
px | 0.743 0.691} 0.640 0.679] 0.743 0.731] 0.746 0.120] 0.677 0.786] 0.724 0.789 
do 0.497 0.458 0.444 ~ 0,049 0.117 0.411 
in — 0.055 0.133 — 0.016 0.412 0.336 0.289 
px | 1.000 0.037] 1.000 0.703] 1.000 0.828} 1.000 0.155] 1.000 0.743] 1,000 0.624 
a |  —-0.761 0.412 0.467 0.405 0.116 ~ 0.015 
16 0.647 0.579 ~ 0.31] ~ 0.901 0.659 0.781 
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Fig. 4. Symmetry shapes for the atomic orbital parts of the occupied symmetry adapted octa- 
hedral site orbitals. | 


it 
As is seen from Table 2 and Fig. 1, the energy bands found by solving the 
different secular equations corresponding to [1], section 3, are lowest in gener: 
for Z=3, but this minimum value varies between Z=2 and 4 for different or- } 
bitals and different wave vectors. For the Z values minimizing the energy the 
interaction range of the corresponding orbitals is, however, sufficiently large to 
give indefinite overlap matrices for many of the secular equations, though boron 
neighbours up to the fourth order of magnitude (24 nearest neighbours) have 
been taken into account. This is also the reason why many of the columns in. 
Table 2 are incomplete. The numerical results obtained, however, give sufficient 
evidence for establishing the existence of minima for the orbital energies for Z 
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values lying between 2 and 4. Far from the minimum points it is found that 
the orbital energies tend asymptotically towards the L approximation. 

The secular equations have been solved in general only for wave vectors with 
high symmetry in the reciprocal space. In the L approximation, however, the 
bands have been given for seven directions of high symmetry including many 
intermediate points (Figs. 2a-e), although the symmetry adapted orbitals are in 
fact somewhat mixed for these points. But, since this mixing enters only in the 
interaction integrals between neighbouring sites, it might be neglected. Fig. 3 
gives a direct comparison between the regions of the energy bands in the L 
approximation and the LM approximation for the best choice of Z. 

The atomic orbital parts of the different occupied symmetry adapted site or- 
bitals are described in Fig. 4. The coefficients for these atomic orbital parts 
occurring in the hybridizations have been given for a few wave vector points 
in Table 3 for the different occupied symmetry adapted orbitals in the LZ and 
LM approximations. 


3. Discussion of the results 


The results given in the previous section show the importance of taking M 
shell orbitals into account. Though the sequences of all the common levels are 
generally equal in the L and LM approximations, the lowest lying orbital energies 
are essentially lowered, when 3d orbitals are taken into account. Strong indica- 
tion was found for the supposition made in [1] that the crystal orbital exponents 
for the M shell electrons should be chosen much larger than for the free atom 
case. It is interesting that the most suitable Z values are so remarkably high 
(Table 2 and Fig. 1). This depends on the ligand field of the crystal. For simi- 
lar reasons a small effect can probably be expected, also when the L shell crystal 
orbital exponents are varied. This fact remains to be investigated. 

The energy bands in the ZL approximation, given by Fig. 2, can be compared 
to the energy bands given by Longuet-Higgins et al. [5], who consider the same 
case. The main difference is that Longuet-Higgins et al. have used several very 
rough approximations and assumptions in addition to those used here. These 
assumptions are as follows: 

1. Two-centre energy integrals are supposed to be proportional to the corre- 
sponding overlap integrals. 

2. 2s and 2p orbital integrals are supposed to be equal. 

3. Only the first order of boron neighbours are considered. 

On account of these approximate assumptions, one cannot expect the present 
results to be essentially in accordance with [5]. The main difference is that some 
of the present bands are conduction bands, while Longuet-Higgins et al. have 
found a complete shell description. Moreover, the sequences of the energy bands 
are somewhat different in the two investigations. 

Fig. 3 is suitable for a comparison between the J approximation and the LM 
approximation. It is seen that the lowest lying levels are lowered, while the 
higher levels are further raised by the introduction of d orbitals. Both the approxi- 
mations have the qualitative results in common, that there are completely oc- 
eupied bonding orbital bands, Itiu, le, and la, and conduction bands, 3t1u, 
lfey, lt2, and 1t,,. In the LM approximation some new conduction bands also 
occur (3t2,, 2t1, and 4e,). These bands are, however, very wide and cannot be 
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given with great accuracy, since they do not exist in the L approximation. The 
main difference from the results of [1] is the fact that there is a small finite 
gap between the fully occupied bands and the conduction bands (see Fig. 3). 

The occupation numbers of the different bands can be roughly estimated by 
considering the number of calculated points of the energy bands with sufficiently 
low energy. A rough description of the configuration can be given by 


(1 tiu)® (1e9)* (1 @1g)? (3 tau)? (1b 29)? (22) 


Thus 12 valence electrons occupy a closed shell and the remaining 8 are con- 
duction electrons moving throughout the whole crystal. A certain very small 
contribution may come, however, from the conduction bands 3teg, 2tiyg, lt and 
lts,. The orbitals 2a1,, 2t:, and 2e, occurring in both the Z and LM approxi- 
mations are antibonding. The results given above form very good verifications 
for the supposition that the crystal is ionized and for the conductivity of the 
bivalent metal borides of type MB,. 

The coefficients of the different atomic orbital parts given in Table 3 indicate 
a remarkably great contribution from the d orbitals, although the corresponding 
numerical values may not be considered as quite reliable, on account of the large 
d overlap. The qualitative results, however, are ascertained and the experimen- 
tally found stability and conductive properties of the CaB, crystal are satis- 
fyingly verified. 
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